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A fluidized granular medium as an instance of the Fluctuation Theorem 
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We study the statistics of the power flux into a collection of inelastic beads maintained in a flu- 
idized steady-state by external mechanical driving. The power shows large fluctuations, including 
frequent large negative fluctuations, about its average value. The relative probabilities of positive 
and negative fluctuations in the power flux are in close accord with the Fluctuation Theorem of 
Gallavotti and Cohen, even at time scales shorter than those required by the theorem. We also 
compare an effective temperature that emerges from this analysis to the kinetic granular tempera- 
ture. 
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Take a fistful of marbles in your hand and shake them 
vigorously. In order to maintain the motions of the mar- 
bles, you will have to continuously supply enough energy 
to them to replenish the energy they lose when they col- 
lide inelastically with each other. This is unlike a state 
of thermal equilibrium, where equal amounts of heat flow 
into and out of the system of particles to the thermal 
bath that sets the temperature of the system. All fluctu- 
ations in the equilibrium state are given by the appropri- 
ate calculation within the canonical ensemble. When this 
equilibrium state is perturbed slightly (for instance, by 
connecting the system between two heat baths at slightly 
different temperatures), the fluctuation-dissipation theo- 
rem specifies the linear response of the system in terms 
of the equilibrium fluctuations. No equivalent framework 
currently exists for large departures from equilibrium. 
However, a recent theorem due to Gallavotti and Cohen 
P, takes an important step in this direction. Inspired 
by an observation made in a simulation of a sheared 
hard-sphere fluid they proved a very general result 
regarding the entropy flux into a system maintained in 
a nonequilibrium steady state by a time-reversible ther- 
mostat. If dynamics in the system are chaotic Q, then 



n(cr T )/n(-ov) = exp(a T T), 



(1) 



where II (ov) is the probability of a fluctuation of am- 
plitude cr T in the rate of entropy production, computed 
over a time r that is much longer than any microscopic 
time-scale of the system. 

The Fluctuation Theorem (FT) embodied in eqn. ||TJ 
above, is a remarkably strong statement, with no ad- 
justable system-dependent parameters. The significance 
and realm of validity of this result are being explored 
theoretically in various directions. The theorem has been 
shown to be true under Langevin dynamics 0, @ , thus 
broadening greatly the range of applicability of systems 
from only Hamiltonian systems. The theorem has been 
shown to be equivalent to the fluctuation-dissipation the- 
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orem in the limit of small forces ■ The FT has also been 
shown 8] for a chain of anharmonic oscillators to produce 
the Green-Kubo formula for thermal conductivity. Re- 
cent work [9( also shows the connection between the FT 
and a new thermodynamic identity [Toj that relates equi- 
librium free energy differences to nonequilibrium work. 
Thus, while eqn. Q may at first sight appear to be a 
rather formal mathematical relationship, it has a very 
nontrivial physical significance. 

Explicit realizations of the FT have been demonstrated 
in simulations of sheared fluids Q , electrical conductiv- 
ity in an array of fixed scatterers shell- models of 
turbulence 12], a model of a granular medium and a 
Burridge-Knopoff type block-spring model 0] . Clear ex- 
perimental demonstrations of the FT have not previously 
been achieved, but the two previous attempts have been 
very instructive in the difficulties involved. Goldburg et 
al. |b| studied the fluctuations in electrical current re- 
quired to maintain in steady-state a liquid-crystal film 
driven into chaotic convectional motion by a constant 
transverse voltage. No direct test of equation was pos- 
sible due to the fact that negative fluctuations in power 
are exceedingly rare in a macroscopic system. Ciliberto 
and Laroche |15| studied temperature fluctuations in a 
Rayleigh-Benard cell arguing that the temperature vari- 
ations at a point in a chaotic flow are representative of 
the heat fluxes to which the FT applies; they side-step 
the difficulty of measuring negative power fluctuations 
encountered in 0] by using as the variable a the lo- 
cal temperature minus the mean temperature in the cell. 
They find that the ratio on the left side of eqn. is 
indeed exponential in the amplitude of the fluctuations, 
but their results are less conclusive on the r dependence. 

Informed by these previous studies, we are led to con- 
sider in this article a granular medium made up of macro- 
scopic spherical beads in which we are able to directly 
measure energy fluxes. Since granular systems typically 
contain much fewer particles than molecular systems, and 
since relevant experimental situations are often strongly 
driven non-equilibrium states, these fluxes naturally un- 
dergo large fluctuations, including the negative fluctua- 
tions that play an important role in the FT. These at- 
tributes of our experimental system allow us to make a 
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FIG. 1: Sketch of the experimental cell. The dashed rectangle 
is a window measuring 10dx21d, fixed in the laboratory frame, 
in which we study the flux of kinetic energy. 



very substantive test of the result in equation (Q. 

We make measurements of power fluctuations in the 
geometry sketched in Figure Q Spherical glass beads 
(diameter d — 1.6 mm; mass m = 5.24 mg) are held 
in a rectangular cage in the vertical plane and agitated 
strongly by vertical vibrations. For the measurements 
we report here, the frequency and amplitude of vibra- 
tion are held fixed at 60 Hz and 2.6d We take video 
frames every 0.5 ms from which the instantaneous po- 
sitions of the particles are determined to a precision of 
0.025c£, and their velocities are thence inferred. We mea- 
sure the time-dependence of the total energy and fluxes 
into a subsystem defined by the window indicated by 
dashed lines in Fig. ^ We quantify the power P(t) sus- 
taining the dynamics in this subsystem by finding the 
sum of two measured quantities. The first is the work 
done on the particles by gravity, which we obtain from 
the difference in potential energy of the particles in the 
window from one frame to the next; this has a zero time- 
average. The second, and much larger, term is the flux 
in kinetic energy through the dashed lines in Fig. [J due 
to particles entering or exiting the window. We neglect 
the contributions to the P{t) from the flux of rotational 
kinetic energy. We note that these data are taken in a 
dilute regime (area fraction ~ 13.8%), so that P(t) is 
well-approximated by measuring only the contribution 
from the streaming term, i.e. ignoring the rare contribu- 
tions to the energy flux from particles at the edges of the 
subsystem that suffer collisions with particles outside the 
subsystem. 

In the subsystem we consider, the fluctuation over a 
time r of the total mechanical energy is AE = D + P, 
where D is the dissipation due to inelastic collisions be- 
tween particles, and P is the power driving the sub- 
system, as described above. In steady state, the time- 
average AE is zero, so that P{t) necessarily has a positive 
time-average value, P, which must exactly balance the 
time-average of the dissipation D. In theoretical discus- 
sions of the FT, the entropy production rate, <r, is iden- 
tified with the phase space contraction rate. Since this is 
not an experimentally accessible quantity, we have chosen 



20 - 



o 

Q_ 

T3 
CD 
N 



-20 - 



4.0 4.1 4.2 4.3 4.4 4.5 
Time [s] 

FIG. 2: A sample (~ 0.6%) of the normalized power trace 
for the subsystem in consideration. The horizontal line shows 
the average power, p = 1. 



to make a correspondence between a(t) and P(t)/T e ff, 
where T e ff is the effective temperature of the system. 
The choice of what temperature to ascribe to the sys- 
tem is not obvious; we can avoid this issue in construct- 
ing the left side of equation 0J, since it is a ratio be- 
tween the probabilities positive and negative values of 
the same variable. However, we must confront this is- 
sue in discussing the right side, since an energy scale will 
be required to make non-dimensional the quantity in the 
exponent. We defer this discussion to later in this article. 

In Figure [3 we display a sample of the time- variation 
of the normalized power, p{t) = P/P to indicate that this 
quantity fluctuates strongly about its mean (p = 1), in- 
cluding making several negative excursions. We empha- 
size that the occurrence of frequent negative fluctuations 
is not because we are operating at small driving forces, 
close to thermal equilibrium; on the contrary, the mean 
kinetic energy is on the order of 10 16 fc^T, where T is 
the ambient temperature that controls the microscopic 
degrees of freedom internal to the particles. 

With p(t) being identified as the variable of interest, 
the version of the Fluctuation relation that we test here 



n(p r )/n(-p T ) = exp(p T TP/T eff ). 



(2) 



Before embarking on the data analysis required to test 
equation J5J, we note a few points that are relevant to 
the implications of any results that emerge from such an 
analysis: 

• The dynamics in the experimental system are not 
time- reversible unlike Q, 0] nor do they satisfy mi- 
croscopic detailed balance unlike the systems of 

Hi- 

• We study an open subsystem of the entire nonequi- 
librium system rather than the global power injec- 
tion. The FT is proven for global fluxes and local 
versions of equation (0 have not been proved ex- 
cept for special cases |l6j | . 
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FIG. 3: Probability distribution, H(p T ) of the binned power 
in time bins r = 0.5,1,2,4,8 and 16 ms. The distributions 
are displaced vertically for clarity. The average power for all 
distributions is p T = 1, indicated by the dashed vertical line. 

• To prefigure the data analysis to follow, we find 
that we are not able to satisfy the condition that 
t is much greater than the dynamical time-scales 
in the system. However, we find that the statistics 
of p(t) are strongly non-gaussian for smaller r and 
therefore provide a very demanding test of equation 
©■ 

In order to construct the left side of equation J5J), we 
first partition the time trace of p(t) into non-overlapping 
bins of duration r and compute the average power 
in these bins, p T (t) = 7 ft p{t')dt' '. The probabil- 
ity distribution Id (p T ) is displayed in figure for r = 
0.5, 1, 2, 4, 8 and 16 ms for a data set constructed from a 
sequence of 170,000 video frames. The data show large, 
non-gaussian deviations from the mean value of p T = 1 , 
with substantial negative tails. The shapes of the distri- 
butions vary with r, thus each distribution provides an 
independent test of eqn. J3J. As noted earlier, even with 
the large data set that we consider, the values of r we 
are able to access are not truly long compared to the dy- 
namical time-scales of the system, as required by the FT. 
For comparison, the mean free time is 4 ms and the time 
for a particle to diffuse across the window is 10 ms. The 
autocorrelation function of p(t) shows weak oscillations 
with a period of 8 ms, which corresponds to the second 
harmonic of the drive frequency: this is due to energy in- 
jection from the top and bottom walls in opposite phases 
of the drive cycle. 

In Fig. QJY we show for the probability distributions of 
Fig. the variation of the ratio ln[H(p T ) /U(—p T )] with 
p T . This ratio is a straight line, as predicted from the 
exponential dependence in eqn. ([Q). The only deviations 
from linearity are the points close to p — on the r = 0.5 
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FIG. 4: (A) Zn[II(p T )/II(— p T )] versus p T P for r ranging 
from 0.5 to 16 ms. (B) ln[H(p T ) Pt)]/t versus p T P 
(P — 356 m 2 s -3 ). The solid line shows the slope of the 
collapsed curves. A dashed line of slope 1/T gra n is drawn for 
comparison. 



and 1 ms data sets. To study the r-dependence implied 
by the FT, in Fig. ^3 we plot the dependence on r of 
the quantity ln[H(p T ) /U(~ p T )]/r. For all values of r, 
we obtain the same slope (similar results were observed 
in the simulations in Ref [lj|)j and for large r, all the 
data sets collapse on a line passing through the origin, 
as predicted by eqn. J5J. This data collapse indicates a 
significant degree of agreement with the FT since each 
of the lines derive from differently shaped H(p T ). Fur- 
thermore the linear dependence of the ordinate in Figure 
4 extends to values of p T well beyond the mean so that 
it seems unlikely that the linearity is merely the leading 
order behavior in a more complicated functional depen- 
dence on p T . 

Finally, we consider the slope of the lines in Fig. 21 
as discussed earlier, had we used the rate of entropy flow 
into the system these lines are predicted to have a uni- 
versal slope of unity. With our choice of variable, viz. 
the power flux, the lines in Fig. 0]have a slope with the 
dimension of an inverse energy/mass, which should be 
given by the inverse of the effective temperature, T e ff, 
of the system. T e ff bears no relation to the ambient 
temperature. A plausible candidate for a temperature 
scale, however, is the so-called "granular temperature" 
Tgran = 1/2 (v 2 ) , a purely kineticall y d efined tempera- 
ture that we have previously studied [2(J, W$ m the same 
experimental geometry. The dotted line in Fig. 0|has a 
slope of 1/T gran : T e ff and T gran are clearly unequal. To 
study how they are related to one other, in Figure [S] we 
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FIG. 5: T e f f (left axis) , defined as the inverse of the slope in a 
graph such as Figure 4B, and T gran (right axis), the granular 
temperature (l/2(i> 2 )), as a function of area fraction. T e // is 
numerically larger but follows a similar trend. 



display the variation of T e / / and T gran as a function of 
the number of particles in the cell, keeping all other pa- 
rameters of the driving mechanism fixed. As the number 
of particles increases, T gran decreases due to an increased 
rate of inelastic collisions. (We make use of the fact that 
the slopes of the lines in Fig. 0]do not depend on r, and 
plot only the values of the slope for r = 1 ms, for which 
more modest statistics are sufficient). As can be seen, 
T e ff and T gran differ in magnitude by about a factor of 
8, however they both decrease similarly as the number of 
particles is varied. 

Prior to the proof of the FT presented in [HQ. an 
even stronger statement was proven by Evans and Sear- 
les in which they showed that eqn. was true at 
all t. This statement, however, applies only to systems 
that start in equilibrium and are perturbed by a time- 



independent force. A recent experiment 011 a c °l~ 
loidal particle pulled through a fluid by an optical trap 
tests this transient fluctuation theorem by examining an 
integral version of eqn. ft|. The transient fluctuation 
theorem of Evans and Searles, however, is not expected to 
apply to the noncquilibrium steady states that we probe 
here [iH . The fact that eqn. appears to hold at 
small r in our experimental situation should be viewed 
as a characteristic of the structure of nonequilibrium sta- 
tionary states that is specific to fluidized granular media. 

As discussed earlier, the experimental system is not 
an idealized instance of the FT: the conditions of time- 
reversibility or detailed balance do not obtain in the ex- 
periment, the variable p(t) is not an entropy generation 
rate, we only study a small, open, subsystem rather than 
the entire system, and we are not able to go to extremely 
long timescales. Despite this, we find excellent agreement 
between our results and the predictions of the FT. This 
gives hope that these ideas can be extended to some sit- 
uations where the conditions of proof are not met. The 
interpretation of the effective temperature in this anal- 
ysis remains an open question [23J. We speculate that 
T e ff is potentially valuable in regimes where the granu- 
lar temperature is no longer expected to be useful. 
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